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Abstract 

This paper develops one of the methods for study of nonlinear Par- 
tial Differential equations. We generalize Sato equation and represent 
the algorithm for construction of some classes of nonlinear Partial Dif- 
ferential Equations (PDE) together with solutions parameterized by 
the set of arbitrary functions. 

1 Introduction 

Nonlinear Partial Differential Equations (PDE) have a great variety of appli- 
cation in both mathematical physics and applied mathematics. The well de- 
veloped tools for study are admitted by the wide class of PDE which is called 
Completely Integrable Nonlinear PDE. These systems reflect many impor- 
tant properties of real physical objects and have been studied in many details. 
So, Korteweg-de Vries equation (KdV) [1, 2, 3, 4], Kadomtsev-Petviashvili 
equation (KP) [4], Camassa-Holm equation [5, 6, 7, 8] and some of its two 
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dimensional generalizations (for instance, [9]), are well known in hydrody- 
namics, where they describe propagation of surface waves on incompress- 
ible fluid without viscosity; Nonlinear Schrodinger equation (NLS) is closely 
used in nonlinear fiber optics, where it describes propagation of electromag- 
netic pulses [10, 11]. Note that some methods for study of the integrable 
PDE do not presume in advance their integrability. For instance, Sato the- 
ory [12, 13, 14], Hirota bilinear method [15, 16, 17, 18], Penleve method 
[19, 20, 21], recent modification of the dressing method based on the alge- 
braic matrix equation [22, 23, 24]. 

In this paper we use the ideas of the ref. [12], where the space of one 
dimensional pseudo-differential (or differential) operators has been used for 
description of the KP hierarchy. We replace the space of one dimensional 
operators with the space of multi-dimensional differential operators, which 
allows one to treat larger class of PDE, in particular, (n + l)-dimensional 
PDE with n > 2. We provide the family of particular solutions, parameter- 
ized by the set of arbitrary functions. Although we have not answered the 
question about complete integrability of the PDE under consideration. In 
comparison with approach considered in the refs. [22, 23, 24], the manifold 
of available solutions is reacher. Trying to make the understanding of this 
article independent on references, we supply all necessary details below. 

General idea is developed in Sec. 2. The useful subspace of n-order multi- 
dimensional differential operators together with related systems of nonlinear 
PDE is considered in the Sec. 3. We represent simple example of (2+1)- 
dimensional equation and generalization of the KP hierarchy in the same 
section. Algorithm for construction of the families of particular solutions is 
discussed in the Sec. 4. Conclusions are given in the Sec. 5. 



2 Space of N- dimensional differential opera- 
tors and related systems of nonlinear PDE 

Hereafter we will use Greek letters to denote multi-index, a = («i . . .chat), 
where N is the dimension of x-space. The dimension of t-space (space of 
additional parameters £ 7 , 7 = (71, . . .7^) ) may be arbitrary. In particular, 
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it may be infinity. By definition, 

N 

x = (x 1 ,...,x N ), \\a\\ = y^Oj, (1) 

i=i 

N Ft Ft 

is the number of all possible values of a with ||o;|| = n. (2) 

All functions are differentiable with respect to their arguments as many times 
as one needs. 

Denote P n the linear space of n-order differential operators of the following 
general form 

n 

W n e P n : W n = w P d ^ ( 3 ) 

where wp are functions of x and t, which will be specified below. Denote Q n 
the number of terms in the operator W n of general form (3): 

n 

Qn = Qn- (4) 

|| Q ||=o 

Space ¥ n is related with Qn-dimensional linear space C^ n of functions ip(x, t) 
with the basis ipi, i = 1, . . . , Q n , satisfying the following condition: 

detail ^0, 0<\\P\\<n, t = l,...,Q n , (5) 

where is Q n x Q n matrix. Indexes a and i enumerate columns and 

rows of this matrix respectively. By definition, dependence on parameters 
x is arbitrary, while parameters £ 7 are introduced in functions ipi by the 
following formulas: 

d 7 ipi = d 1 ^. (6) 

Suggestion 1. Let parameter Q n be defined by the eq.(4). For any a 
such that ||a|| — n and any space C Qn_1 there is set of Q n -i uniquely defined 
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functions Wp such that space C^ n_1 forms Q n _i-dimensional linear subspace 
of solutions for the following linear differential equation: 

n—1 

S a ^j = 0, S a = d a + w a , w a = W P P e pn - ( 7 ) 

11011=0 

A Let us fix parameter a with | \a\ | — n and basis ipi, i — 1, . . . , Q n -\ of 
space C^" -1 . One can write the formal system of differential equations (7) 

n—1 

S a if>i = d a A + w P^i = °> i = 1, • • • , Qn-i- (8) 



Let us consider the system (8) as nonhomogeneous system of algebraic equa- 
tions on the functions w*. The matrix of this system, (0 < \\/3\\ < 
n — l,i = l,..., Q n -i), is nondegenerate due to the condition (5). So that 
this system has unique solution, i.e. functions Wp are defined uniquely in 
terms of functions ipi and their derivatives. Owing to this, Wp are functions 
of both sets of variables, x and t. By construction, we have determined coef- 
ficients Wg in the system of differential equations (7), for which given set of 
functions ipi (i — 1, . . . , Q n -i) forms basis for certain solution subspace. ▲ 

Equations of the form (7) with all possible values of parameter a such that 
||q;|| = n, form system of differential equations. One can show that its 
coefficients w a satisfy some system of nonlinear PDE. For this purpose, one 
needs to fix multi-index 7, ||7|| = n 7 , differentiate the eqs.(7) with respect 
to i 7 and study the result: 

^ + s a er^i, = o, || 7 || = n 7 . (9) 

Suggestion 2. For each particular pair of parameters a (\\a\\ = n) and 
7 there is set of Q% differential operators B iae = Y^\\s\\=q v i^sd 5 G P™ 7 with 
coefficients expressed in terms of functions Wp and their derivatives such that 
operators M ia 

M ia = + S a ff> + B ^ ( 10 ) 
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belong to the space P n , i.e. M ia can be written in the form 



ra-l 



M ja = a T d ^ (11) 



where a 7 / 1 are some functions of uj^ and their derivatives. 

A The first term in the operator M 7a , = ~Q~ ' belongs to 

7 ll/3||=0 7 

the space P n_1 . Since .S" G P n , the second term in the eq.(10) belongs to 
the space P n+r V Let B iae = £!|?||=o u 7arfd* G P n ^ be n 7 -order differential 
operators with arbitrary coefficients. Then one can write 



B iae S e = b, a ^ + R n -\ fi^Gf- 1 , (12) 

\\e\\=n \\fl\\=n 

b-ya/3 = V 7 ae(/3-e) + (13) 

/(w |, \\S\\ > H/3 — e\ |; coefficients of operators W a ), 



\\(3\ \ —n,...,n + n. 



7- 



This means that J]||e||=n B ia£ S e G P n+n T and has all arbitrary coefficients 
ahead of the derivatives of the order grater or equal to n. Let us choose 
these coefficients in such a way that 



S a d~<+ B iaeSA GP"- 1 . (14) 



which provides condition M~ /a G P n_1 . Due to the relation (14), coefficients 
v iae s are expressed in terms of the coefficients of the operators W a , \ \a\ \ = n. 
A 

By construction, operators M ia generate the following system of Q n -i x 
Q" linear differential equations: 

n— 1 

M ia ^i= ^ aj x d (3 ^ i = 0, ||a||=7i, i = 1, . . . , Q n -\- (15) 
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Let us regard the system (15) as linear homogeneous system of algebraic 
equations on coefficients a J* with nondegenerate matrix (0 < < 

n — 1) due to the condition (5). This system has only trivial solution, i.e. 



From another point of view, coefficients d^ a have been defined in terms of 
coefficients of the operators W a and their derivatives, i.e. the system (16) is 
system of nonlinear PDE on functions u/g, generated by the space of (n — 1)- 
order differential operators W a with parameter t 7 introduced by the eq.(6). 
The complete set of systems of nonlinear PDE, related with all possible 
parameters t 7 , introduced by the eq.(6), forms hierarchy of PDE associated 
with given space P n_1 . 

We call the equation (15) written in the form 



generalization of Sato equation. The system of equations (8) serves for con- 
struction of particular solutions to the system of nonlinear PDE (16), which 
will be used in Sec. 4. 

The system (16) consists of Q n -i x Qn equations. Equations of the system 
(16) with any particular 7 recursively relate all Q n -i x Qn coefficients 
(\\a\\ = n, \\(3\\ = 0, 1, . . . , n — 1) of the operators W a and compose the 
complete system of equations for whole set of these coefficients. Disadvantage 
of this system is that its structure depends on n. In the next section we 
consider the subspace of the space P n_1 , which allows one to generate the 
systems of PDE which do not depend on the value of parameter n, if only 
this value is big enough. These systems of PDE admit the class of solutions, 
which depends on the set of arbitrary functions of single variable. 



a 



0. 



(16) 




(17) 



3 About reductions 



Let V\ 



in 

U2,...,n N 



CP",n = J2i=i n ii b e defined as follows: 



n 




(18) 



0ll=o 



(3 N ), f3=((n-(3 1 -l)p 2 ...p N ) 
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i.e. parameters rii, i — 1, . . . , N, mean the order of the highest derivative with 
respect to x« in the differential operators W n . Emphasize that we don't put 
any restriction on the order of the derivative with respect to x 1 , i.e its highest 
possible order in operators W n equals n. Index with hat is introduced for 
convenience of representation of system of nonlinear PDE. All multi-indexes 
(without hat) in this section have the following structure: 

a = {oi\ ■ ■ ■ oln) with cti < rii for i > 1. (19) 

Let us introduce the set of equations, which specifies the dependence on 
Xi, i > 1: 

^ + V = «9 a<l V, a® = (a®... a <$), a? = 0, (20) 

where tp G C® n . Recall that Q n is the maximum number of terms in operators 
W n G P" 2 ... n . Due to the relations (20), one can use variables (i > 1) 
along with t 7 for construction of the nonlinear eqs. (16). One needs just 
replace <9 7 with <9j in formulas (9-17). This fact will be used in sections 3.1.1 
and 3.1.2. 

3.1 Examples of nonlinear PDE related with space Vf 

In this section we will use two-dimensional x-space, x = (xi,x 2 ), Wj G V^ 1 , 
j = 1,2,.... Maximum order of derivative with respect to x 2 in the operators 
Wj equals one: n 2 = 1. Because of that, the set of operators S a consists of 
two operators Si = S n0 and S 2 = S( n -i)i f° r an Y fixed parameter n: 

n—2 n—l 

S 1 = d^ + W 1 , Wi=J2 u [n - k - 2)1 d$d 2 + M(„- fc -i)o^, (21) 

fc=0 k=0 
n—2 n—l 

S 2 = d r r 1 d 2 + W 2 , W 2 = Y,nn-k-2)id ! {d 2 + J2v(n-k-i)od k 1 . (22) 

k=0 k=0 

These operators introduce two differential equations (7) on the function -0: 

= 0, S 2 il) = 0. (23) 
It is simple to define the number of terms Q n -i in the operators Wf 

Q n _i = 2ra-1. (24) 
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Next, one needs to fix (2n — l)-dimensional basis ipi (i — 1, . . . , 2n — 1) in 
subspace C 2n_1 and write formally the system (8) 

S^A = 0, (25) 
S 2 rl>i = 0, i = l,...,2ra-l, (26) 
Condition (5) now reads 

_ _ / I' (27) 

Each of the systems (25) and (26) should be treated as system of algebraic 
equations on coefficients and % respectively. Due to the condition (27) 
these equations have unique nontrivial solution. 

3.1.1 Example 1: (2+1) dimensional PDE 

In this example we take the eq.(20) in the following form 

%1> = (28) 

and use notations 

x = xi, ti = x 2 , t 2 = t a , a = (20). (29) 

Function ip depends on variable t 2 due to the formula (see eq. (6)) 

d t2 i> = dty. (30) 

Equations (17) related with parameters t\ = x 2 and t 2 have the following 
form: 

dS 

Mtj = tt^ + Sjd, + BijiSi + B lj2 S 2 = 0, j = 1, 2, (31) 
ox 2 

dS 

M 2j = 1 - L + Sjdf + B^S, + B 2j2 S 2 = 0, j = 1,2, (32) 
or 2 

where B^ are differential operators of the next form 

#iii = voo, £n2 = (-uoo + voi) (33) 

B 121 = -1, .B122 = -vqo, (34) 

B 2U = 2u 00x - d\, B 212 = 2u 01x , (35) 

-B221 = 2v 00x , B 222 = 2v 01x - df, (36) 
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which provide the following structure of the operators (see (11)): 

n— 2 n—1 

My = X>?iW + £ a 2,af' i,j = l,2 (37) 

s=0 s=0 

and 

4 = (38) 

(compare with derivation of the eq.(16)). 

Below are several equations from the list (38) which compose the complete 
system of equations: 





-1)0 


= «oi + v 00 v i - via + u 00x2 - 


- ^00* = 0, 


(39) 


n 

a (n . 


-2)1 


= U 10 + U 01 V 00 - U 00 V 01 + Vqi 


- v n + u 01x2 - v 01x = 0, 


(40) 


all 


-1)0 


= -u 00 - VqO + ^01 + V 00 x 2 = 


o, 


(41) 


all 


-2)1 


= -U 1 - VqoVqi + V 10 + V 0ix2 


= 0, 


(42) 


all 


-2)0 


= —Uio - VooV W + Vu + v Wx , 2 


= 0, 


(43) 


afl 


-1)0 


= u 00t2 + 2u ou 00x + 2v 00 u i x 


— 2-uioz ~ uqo xx = 0, 


(44) 


< 


-2)1 


= u 01t2 + 2uoiu 00x + 2v iu 01x 


— 2uu x ~ u 01xx — 0, 


(45) 


< 


-1)0 


= V00t 2 + 2m 00^00x + 2^00^01^ 


— 2v 10a . — v 00xx = 0, 


(46) 


afl 


-2)1 


= v 01t2 + 2m oi^oox + 2v iVoix 


- 2v llx ~ VQlxx = 0, 


(47) 



We leave without proof the statement that this system remains the same for 
any value of parameter n, which happens due to the special structure of the 
subspace V^ 1 - 

To simplify the system (39-47), let us solve eqs. (41,42,43) with respect 
to Moo, Wen an d iiio respectively, substitute the result in the rest of equations. 
Then eqs. (46) and (39,40) result in the following system of equations on the 
functions u = v 00 , v = v ±, and w — v w : 

u t2 — 2u 2 u x + 2vu x + 2u X2 u x + 2uv x — 2w x — u xx = 0, (48) 

2uu X2 ~\~ 2v X2 -\- u X2X2 u x 0, 
-2vu X2 + 2w X2 + v X2X2 -v x = 0. 

Simple example of particular solution for this system will be constructed in 
the Sec. 4. 
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3.1.2 Example 2: (3+l)-dimensional generalization of KP hierar- 
chy 

In this example we use the following notations: 

x = x u h = X2, ti = t a(i) , z = 2,3, a {2) = (11), «( 3 ) = (30). (49) 
Dependence on ti — X2, t 2 and £3 is given by the following equations: 

8^ = 8^, 8^ = 8^, 8^ = 8^, (50) 
Equations (17) have the following form: 
8S- 

M y = 7: - + Sjtit + ByiSi + B lj2 S 2 = 0, 3 = 1, 2, (51) 
0x2 

(99 

M 2j = + + S 2jl S' 1 + B 2j 2S 2 = 0, j = 1, 2, (52) 
<99- 

M v = + S A + BzjiSi + %2^2 = 0, 3 = 1, 2, (53) 
One has the following expressions for differential operators B^: 

Bm = (-M01+U00), (54) 

£112 = (-M00 + ^01) - <9i, (55) 

£121 = («oo - voi) - di, (56) 

#122 = (uoi-voo), (57) 

-B211 = (uoouoi - uu - v 00 voi + ^10 + Zvoox) + (-U01 + v 00 )d 1 , (58) 

#212 = («oi - ^10 - ^01^00 + ^00^01 - «oi + U H + 2v 01o;) + (59) 

(-w 00 + ^oiH - <9i, 

-S221 = (-«oo + u w - U iV 00 + ^00 + ^oo^oi - Vll + 2u 00x ) + (60) 

(^00 - v 01 )di - 8f, 
B222 = (-(u 00 u i)+uii + v 00 v i-v w + 2u i x ) + (u i-v 00 )8i, (61) 

-B311 = -3(1x00^000; + ^oo^oia; - ui 0x - u 00xx ) + 3u 00x 8i - df, (62) 

#312 = -H u oiUoox + VoiUoi x -u 1 i x -u i xx ) +3u i x 8i, (63) 

-S321 = -3(u 00 voo x + v 00 v i x - vi 0x - v 00xx ) + 3v 00x 8i, (64) 

S322 = -3(^01^00^ + foi^oix - v Ux - v oixx) + 3v i x di - df. (65) 
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With given operators Bijk, equations (51-53) take the form (37) with % = 
1,2,3 and j = 1,2. Write down several nonlinear equations, generated by 
each of the eqs. (51-53). 
From the eq.(51): 

-(uooUoi) + U U + VqoVqi - V 10 + UoOx 2 - V 00x = 0, (66) 
-«0! + UlO + MOI^OO - ^00^01 + ^01 ~ V U + ( 67 ) 
u 01x 2 ~~ V 01rr — 0, 

From the eq.(52): 

a (n-l)0 = W 00^01 - ^01^10 - ^00^11 + «21 + ^01^00 ~ ^01^00 _ ( 68 ) 

voo^oi + voi^io + voo^n - ^20 + «oot 2 - u m u mx + 

^00%0:r + W00^00:r + V lV00:r + 2v 0^01x - 
a (n-2)l = "00^01 - 2w lMll + «20 + Mll^OO + ^01^01 ~ (69) 

Mio^oi - 2m 01 w oWoi + Moo^oi - v l\ + M oi^io - 

WOO^ll + 2^01^11 - v 2 i + u 01t2 - U iU i x + V 00 U i x + 

2u 01 v 0Qx - u 0Q v 01x + 3v 01 v 01x - %v llx - v 01xx = 0, 



From the eq.(53): 

Uoot 3 - 3«oo w oox + 3u 10 u 00x - 3u 0l v 00 u 00x + 3u 00 2 x - (70) 
3u 00 v 00 u i x - 3v 00 v iu 01x + 3v w u 01x + 3u 00 u Wx + 

3VqoUh x — 3U20x + 3Uqi x V()0x + 3UqqUqq xx + 
3VqqUqi xx — 3Uiq xx — Uqq xxx = 0, 



Let us show that the above lists of nonlinear PDE can be reduced to 
the integrated KP. First of all note that the reduction d 2 = d 1 leads to the 
following identities: 

Vij = Uij, un = Uio, for all % and j. (71) 



n 11 

U (n-1)0 
a 11 

a (n-2)l 



a 



31 

(n-l)O 
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Then equations (66) and (67) becomes identities, while the system (68) - (70) 
results in three equations ( u = woo, v — u w, w — u 2o) 

u tl + 4uu x - 2v x - u xx = 0, v tl + 4vu x - 2w x - v xx = 0, (72) 
Ut 3 — u xxx — 3v xx — 3w x + 6ul + 6uu xx + 6vu x + 6uv x — 12u 2 u x = 0. (73) 

Eq.(73) results in the integrated KP after elimination of v and w: 

1 3 

u t3 - -^Uxxx - -d~ l u hh + 3u 2 x = (74) 

For this reason we call the above system (66) - (70) (together with the whole 
list of nonlinear PDE related with different variables t 7 ) the generalization 
of KP hierarchy. 



4 Construction of particular solutions 

For construction of particular solutions to the nonlinear PDE which appear 
in this paper, one needs to choose the functions fy, i = 1, . . . , Q n -i which 
satisfy the conditions (5) and (6). One of the following Fourier integrals is 
suitable for this purpose. 

For the equations from the Sec. 2: 



= / Cj(k) exp | kiXi + ^^u; 7 (k)t 7 J dk, (75) 

J \i=l 7 / 



(ki, . . . , k N ), i — 1, . . . , Qn-i- 



Here uj 1 depends on k due to the dispersion equations associated with the 
eqs.(6) 

u-y — k 1 . (76) 

In this case ipi are arbitrary functions of all variables x«. 
For the equations from the Sec. 3: 



ipi = J Cj(fci) exp ( kiXi + ki{ki)Xi + o> 7 (fci)t 7 j dfa, (77) 



i=2 7 

1 1 , . . . , Q n ~l- 



12 



Functions u 7 (ki) and ki{k\), i > 1 represent the dispersion relations associ- 
ated with eqs.(6) and (20): 

N (0 

uj, = k\ kr +1 = ]Jkf , af } = 0, i = 2,3,...,N. (78) 

3=1 

Here tpi are arbitrary functions of single variable x\\ Cj(k) and Ci(k\) in for- 
mulas (75) and (77) are such functions of argument (s) that condition (5) is 
satisfied. 

After this, one needs to solve the system of algebraic equations (8) for 
coefficients which are the solutions of appropriate system of nonlinear 
PDE (16). 

For instance, let us construct the particular solution for the system (48). 
In this case the equation (77) should be written in the form 

if>i — J Ci(k 2 ) exp {k\xi + k 2 x 2 + k 2 t ± ) dk 2 , i — 1, . . . , Q n -i- (79) 

To construct the simple solution, let Q n -i = 3 (or n = 2 due to the eq.(24)) 
and take the following expressions for i/jf 

fpi = cu + c 2i exp(kiXi + kiX 2 + kft±), % — 1, 2, 3. (80) 

Then the systems (25) and (26) becomes 3x3 algebraic matrix equations 
for the coefficients and v ij respectively, which can be solved directly. One 
can see that the condition (27) is satisfied for our choice of functions ipi. We 
give only expression for the function u = Voq: 

s 1 E 1 E 2 + s 2 E 1 E 3 + s 3 E 2 E 3 + s 4 E 1 E 2 E 3 

u = , (81) 

p±E 1 E 2 + p 2 E 1 E 3 + p 3 E 2 E 3 + p 4 E x E 2 E 3 

Ei = exp(k^x + + kftx), 

si = c 13 c 21 c 22 kik 2 (k 2 2 - kf), s 2 = c 12 c 21 c 23 k 4 k 3 (kl - kj), 

S3 — c \\ c 22C-2 , ik 2 k 3 {kj > — k 2 ), 

s 4 = c 21 c 22 c 23 k 1 k 2 (k 2 -k 1 )(k 1 -k 3 )(k 2 -k 3 )(k 1 + k 2 + k 3 ), 

Pi = cisCjjiczjfciAfe^i - k 2 ), p 2 = ci 2 C2iC23fci*;3(A;3- fci), 

P3 = c u c 22 c 23 k 2 k 3 (k 2 - k 3 ), p 4 = c 2 ic 22 c 23 {k 1 -k 2 ){k 1 -k 3 )(k 2 -k 3 ), 
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If all coefficients pi in the denominator of this solution have the same sign 
(this can be arranged by the appropriate choice of parameters q and fcj), 
then the solution has no singularities for all values of independent variables 
and describes inelastic three-kink interaction. 

5 Conclusions 

We represent an example of multidimensional generalization of Sato equation 
(17), written in terms of pure differential operators. We give algorithm for 
construction of different multidimensional hierarchies which admit the special 
family of particular solutions. Namely, solutions are parameterized by the 
set of arbitrary functions of either one or several variables X{. Different 
subspaces of general space (3) result in specific systems of nonlinear PDE. 
Two examples of this kind are given in the Sec. 3. 

The represented approach might be considered as the mixture of Sato 
theory and Lax representation. In fact, each equation of the system (17) is 
nothing but the compatibility condition of the overdetermined linear systems 
(6) and (7). For each particular 7, the equation (17) gives the complete 
system of nonlinear PDE (16) on the coefficients of the operators W a . 
But the structure of this system depends on n (order of differential operators 
S a ) in general case. The examples of systems of nonlinear PDE which do not 
depend on n are given in the Sec. 3. 

Author thanks Professor A.H.Zimerman (IFT UNESP, Sao Paulo, Brazil) 
for initiation of this work. The work is supported by RFBR grants 01-01- 
00929 and 00-15-96007. 
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